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Abstract 

We study higher order approximations in the renormalization group 
approach to matrix models. We use constraint equations on the free 
energy resulting from a freedom of field redefinitions and obtain the ef- 
fective beta function for a single coupling constant to the fifth order. The 
fixed point and the string susceptibility exponent are shown to approach 
the values obtained in the exact solution as the order of approximations 
becomes higher. 
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Matrix models [1] give a discrete version of two-dimensional quantum gravity. 
By studying the double scaling limit of the matrix models [2], one can obtain the 
exact solutions of two-dimensional quantum gravity coupled to minimal conformal 
field theories with central charges c < 1. These results are consistent with the 
results of a continuum version of two-dimensional quantum gravity [3]. For the 
central charge c > 1, we can construct matrix models [4], but cannot obtain exact 
solutions. Moreover, for those cases we do not know whether there is a continuum 
theory or not. Brezin and Zinn- Justin [5] have proposed a renormalization group 
approach to the matrix models in order to understand c > 1 cases. One can check 
the validity of this approach by applying it to c < 1 theories, since the exact solutions 
are known for these cases. In ref. [5], they obtained reasonable results for the first 
order approximation in the case of c < 1. The renormalization group approach has 
been further studied in several papers [6-10]. 

The purpose of this paper is to study higher order approximations in the renor- 
malization group approach to the matrix models. We consider a one-matrix model 
corresponding to the case c = 0. It was noticed in ref. [9] that one has to use a free- 
dom of field redefinitions to obtain higher order results. The exact beta functions of 
the vector model were obtained in this way [9]. Although it is difficult to obtain ex- 
act results for the matrix models using the renormalization group approach, we can 
compute beta functions perturbatively in a coupling constant. We derive constraint 
equations on the free energy resulting from the freedom of field redefinitions. Using 
these equations we obtain the effective beta function for a single coupling constant 
to the fifth order. The fixed point of the beta function and the string susceptibility 
exponent are shown to approach the values obtained in the exact solution of the 
model as the order of approximations becomes higher. 

We shall study the following one-matrix model: 
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where <pN is &n N x N hermitian matrix. The terms go^mtr^ (m > 1) generate the 
usual discretizations of random surfaces and the terms gk,itr<j)%tr<j) l N (k > 1,1 > k) 
generate two random surfaces touching each other at one point [11,12]. This model 
was proposed in ref. [11] to include effects of higher order intrinsic curvature terms. 
Following ref. [5], by integrating out one row and one column of (N + 1) x (TV + 1) 
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matrix (pN+i, we will see that the partition function Z^{g) fulfills the equation 



ZN+i(g) = T N (g)Z N (g + Sg), 



(2) 



where T^ig) is a factor which is a function of N and g. From this relation we have 
a differential equation 
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Now, it is convenient to introduce rescaled coupling constants [9], 
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By this rescaling, the quadratic term in the potential becomes the standard form 
\trcf> 2 N , i.e., (?o,2 = 1- It is clear that we have a relation 
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where Z^{g) and Z^{g) are shorthand notations of Zjv(<7fc,z {k > 1, I > k), (?o,m (m > 
1)) and Z N (g k;i (k > 1, / > fc), £0,2 = 1, #o,m (m = 1,3, • • •)) respectively. We de- 
note the partial derivatives with respect to g~k,i (k > 1,1 > k), go,m (m — 1,3, • • •) 
and <7o,2 by \g, and those with g^i (k > 1, I > k) and (? ,m (m > 1) by \ g . They are 
related as 
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If we define the free energy 



F(N,g) = -j^lnZ N (g), 



the renormalization group equation (H) becomes 
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Now we shall consider the simplest case of a cubic potential 
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5o,2^0, 5o,3^0, 5o, m = (m = l,4,---), g k ,i = (k>l,l>k). (14) 

The shifts 5g in eq. @ can be obtained exactly as follows. The partition function 
Ziy+i{g) is defined by 
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The matrix 4>n+i is parametrized in terms of a sub-matrix 0tv, a complex vector vn 
and a real number a: 
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Since all the terms involving a are of relative order 1/N [5], we can set a = 0. Hence 
we have 
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Since the ujv and v* N integrals are Gaussian, we can evaluate them exactly. After 
rescaling of 4>n we obtain 
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where 



<^*M = (A; > 1, I > k), 5g 0t2 = |#o,2 - 
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The factor in eq. (Q) is given by T N (g) = { — (N + 1)^0,2}^- We can also consider 
cases of other potentials, but then the vn and integrals cannot be evaluated 
exactly. From eqs. (JTTJ) and ([19D, we have the following renormalization group 
equation for the cubic potential 
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where /3 functions are 
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We find that even if one starts from the cubic potential, after integrating one 
low and one column one obtains infinite number of other higher order terms trcf)™. 
To do higher order calculations, we must take into account these higher order terms. 
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However, one may be able to put the potential back to a cubic form by appropriate 
redefinitions of 4>n- Such a freedom of field redefinitions was used in ref. [9] to 
obtain the exact beta functions in the vector model. 

The freedom of field redefinitions can be expressed as constraint equations on the 
partition function. These constraints were discussed in refs. [13, 14] and were shown 
to have a Virasoro like structure. In particular they were discussed before taking 
a continuum limit in ref. [13]. They can be obtained by a change of integration 
variables 

4>n — > 4>n + S(j) N (23) 
in the partition function. The transformation considered in ref. [13] are 
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where e is an infinitesimal parameter. We call the constraints corresponding to the 
change of variables ([24]) as L n (n > — 1). In the case of our model ([]]) the constraints 
for the free energy F(N, g) are 
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L n (n > 3) conditions: 
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Here we have set all the coupling constants to be zero except 50,2 and g ,3 as i n ec L- 
fll4f) after the differentiations. 



To do higher order calculations we will need constraints derived from more gen- 
eral transformations of 0jv- Let us consider changes of variables 



6<I>n = jj<P n N +1 tr<f N (n = -1, 0, 1, p > 1). 
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The constraints L niP (n — — 1, 0, 1, p > 1) corresponding to these changes of 
variables are found to be 
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We could consider more general transformations of the form fl30| ) with n > 2, p > 



1 or 



tr^ tr<f% tr<ffi 
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However, these transformations induce terms tr(j) kl tr(j) k2 tr(j) k3 • • • (1 < &i < k 2 < 
k 3 < ■ • •) in the potential. So we have to begin with a more general potential 
than ([I]), which includes terms of the form tr<j) kl tr(p k2 tr(j) k3 • - -. To the order we will 
consider in this paper we do not need to use such general transformations. 

Let us discuss implications of the above constraint equations. The L_i condition 
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F. This means that the term 



)n in the potential can be eliminated by an appropriate field redefinition. The 
meaning of the L condition (EBT) is clear from the definition of the free energy 
F(N,g). The L\ condition (27) shows that 
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using L n>p conditions 0), these nonlinear terms can be rewritten by 

and terms higher order in ^o,3- 

We can now calculate the fixed point <7g 3 and the string susceptibility exponent 
7x using the above results. We assume that g 0j3 is small and use a perturbation 
theory in g 3 . This assumption is justified if a value of the fixed point obtained in 
this way is small. 

We begin with the first order calculation, which taking into account up to and 
including terms of order c/q 3 in the renormalization group equation fl2"Dp. Since 
A),m { m > 4) are of order g™ 3 , we can ignore them in the renormalization group 
equation. Using the L_i condition, we can rewrite /?o,i term as the same form as 
/5o,3 term. Thus we obtain an effective renormalization group equation with a single 
beta function (3^/^ 
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The fixed point #0,3 is determined by ffij (5^3) = and the string susceptibility 

exponent is given by 71 = g, 2 . as in ref. [5]. From eq. (|37|), we find the fixed 
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and the string susceptibility exponent 
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On the other hand, from the exact solution [2] we know that the double scaling limit 
is achieved near the critical value 
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with the string susceptibility exponent 
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This first order calculation is essentially equivalent to the first order calculation in 
ref. [5] for the matrix model with quartic potential. 

Next, for the second order approximation we take into account up to and includ- 
ing terms of order g^ 3 . Ignoring /3o, m (m > 5) and using the L_i and L\ conditions 
we obtain the effective renormalization group equation 
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Here we have ignored a term of order g^ 3 . From eq. fl43|) , we find the fixed point 
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and the exponent 
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(45) 



At the third order approximation we take into account up to and including terms 
of order g^ 3 in the renormalization group equation (ppp. By using the L_i, L 1; L 2 , 
and L_ 12 conditions we have the effective renormalization group equation 
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Although a nonlinear effective beta function j3 2 ^ 2 ^ appears in eq. (Pf), we may 
neglect it because it is of order g^ 3 . Then from eq. (|47D we find the fixed point 
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and the exponent 
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Similarly we can continue to improve the approximation. At the fourth order, 
using the L_ 1; Li, L 2 , L 3 , L_ 12 , L 2 and L 3 conditions and ignoring terms 

of order gl 3 we obtain the effective renormalization group equation 
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At the fifth order, using the L_i, L\, L 2 , L 3 , L 4 , L_ 12 , £-1,3, L Q 2 and L ,3 

conditions and ignoring terms of order g^ 3 , we obtain the effective renormalization 
group equation 
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Thus, as the order of approximations becomes higher, resulting values of the fixed 
point and the string susceptibility exponent approach the exact values (|40D , fl4ip . To 
carry out much higher order calculations, we have to use constraints corresponding 
to more general transformations (|30|) with (n > 2, p > 1) and (|3~5|) starting from 
more general potential than ([[]). 
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